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ABSTRACT 

When extracting the weak lensing shear signal, one may employ either locally normalized or globally nor- 
malized shear estimators. The former is the standard approach when estimating cluster masses, while the latter 
is the more common method among peak finding efforts. While both approaches have identical signal-to-noise 
in the weak lensing limit, it is possible that higher order corrections or systematics considerations make one 
estimator preferable over the other. In this paper, we consider the efficacy of both estimators within the context 
of stacked weak lensing mass estimation in the Dark Energy Survey (DES). We find the two estimators have 
nearly identical statistical precision, even after including higher order corrections, but that these corrections 
must be incorporated into the analysis to avoid observationally relevant biases in the recovered masses. We 
also demonstrate that finite bin-width effects may be significant if not properly accounted for, and that the two 
estimators exhibit different systematics, particularly with respect to contamination of the source catalog by 
foreground galaxies. Thus, the two estimators may be employed as a systematics cross-check of each other. 
Stacked weak lensing in the DES should allow for the mean mass of galaxy clusters to be calibrated to « 2% 
precision (statistical only), which can improve the figure of merit of the DES cluste r abundance experimen t 
by a factor of ~ 3 relative to the self-calibration expectation. A companion paper dSchmidt & Rozol 120101) 
investigates how the two types of estimators considered here impact weak lensing peak finding efforts. 
Subject headings: cosmology: clusters 



1. INTRODUCTION 

Upcoming large photometric surveys such as the Dark En- 
ergy Survey (DES jj, Pan-Starrfl and the Large Synoptic Sur- 
vey Telescope (LSST£| will find hundreds of thousands of 
clusters over large fractions of the sky. These samples hold 
the potential to be our most pow erful tool in understanding 
dark energy dAlbrechtet al. 2006), and appear to be a neces- 
sary component of any set of observables that wishes to distin- 
guish between dark energy and modifie d gravity approache s 
for explaining an accelerating universe (Sha piro et alj|2010h . 
Indeed, even today cluster abundances provide some of 
the tightest constraints on t he amplitude of the low red- 
shift matter power spectrum dMantz et al. 2008; Henry et ail 
120091: IVikhlinin et aLll2009t iRozo et al.ll2010h. and inter esting 
constrains on mod i fications to gravity dSchmidt et al.l l2009t 
iRapetti et al]l20Tot ILombriser et alJl2010h . Not surprisingly 
then, realizing the promise of galaxy clusters as a cosmolog- 
ical probe is of paramount importance for understanding the 
physics driving the current phase of accelerated expansion of 
the universe. 

The most important obstacle that cosmological applications 
of cluster surveys must overcome is the calibration of mass- 
observable relations. That is, the Cold Dark Matter (CDM) 
paradigm of structure formation allows us to predict the abun- 
dance of galaxy clusters as a function of mass, whereas em- 
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pirically we can only recover the abundance of galaxy clus- 
ters as a function of some observable X that correlates with 
mass. Consequently, cosmological investigations of cluster 
abundances require that we carefully calibrate the probability 
P(X\M) that a halo of mass M is included in a survey as a 
cluster with observable X. This is a problem that is particu- 
larly difficult for photometrically selected cluster samples, as 
there is little theoretical understanding of the relation between 
a cluster's galaxy content and its total mass. 

Self-calibration is an elegant attempt to overcome this 
difficulty. In this approach, one simply parameterizes 
P(X\M), and fits for the corresponding parameters relying on 
cluste ring information dLima & Hul 12004 120051: IHu & Cohnl 
120061). and/or the evolut i on of abundances w ith redshift 
dMajumdar & Mohri 12001 iGladders et alj 120071) . One may 
further enhance such sel f-calibration techniques by relying 
on multiple mass tracers dCunhall2009l) . and these approaches 
are expected to be very successful in improving our under- 
standing of dark en ergy relative to other dark energy probes 
dCunha et al.ll2009h . Nevertheless, it is expected that careful 
a priori calibration of the mass-observable relation of galaxy 
clusters from targeted follow-up observations can further en- 
hance the utility of cluster samples ov er and above what 
can be achieved through se lf-calibration dMajumdar & Mohri 
l27Ml200l:rWu et al.ll2010h . 

One way of empirically calibratin g cluster masses is 
through cluster weak lensing stacking dSealfon et al] 120061: 
Johnston et al.l2007atlSheldon et alj2009l:lMandelbaum et all 
20081: iLeauthaud et al.ll2010t IWhite etaLlhoiOh . This tech- 
nique relies on coherently adding the weak lensing signal of 
galaxy clusters at fixed observable in order to estimate the 
mean mass of the stacked galaxy clusters. Relative to esti- 
mating the weak lensing mass of individual clusters, cluster 
stacking has the significant advantage of allowing us to detect 
the lensing signal at significantly lower cluster masses than 
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would otherwise not be possible. Moreover, averaging over 
many halos dramatically reduces the impact of weak lens- 
ing projection effects due to non-correlated structures along 
the line of sight. While such measurements do not in any 
way constrain the scatter of the observable-mass relation, the 
effectiveness of cluster stacking on improving cosmological 
constraints from photometrically selected cluster samples is 
evidenced by the dramatic improvement that this measure- 
ment produced on the cosmologica l constraints derived fro m 
the SDSS maxBCG cluster sample dRozo et al.H2007L 120101) . 

Here, we forecast the precision with which cluster mass 
and concentration can be measured in a DES-like photomet- 
ric survey, and we discuss a variety of systematics associated 
with this measurement. Specifically, motivated by the large 
lensing bias signa l expected in the wea k lensing power spec- 
trum in the DES (ISchmidt et al]l2009*al) . we revisit the ques- 
tion of whether lensing bias ca n have a significant impa ct on 
weak lensing mass calibration (M andelbaum et aDl2005h . We 
also consider finite bin-width correction to density contrast 
estimates, and finally, we discuss how photometric redshift 
source selection can impact the recovered weak lensing sig- 
nal. 

This paper is one of two companion papers; the second pa- 
per considers the problem of identifying weak lensing peaks 
in large photometric surveys and the impact of lensing bias 
on weak lensing peak finding. An interesting by-product of 
having performed a simultaneous study of these two distinct 
subjects was the realization that in attacking these problems, 
different approaches for extracting the weak lensing signal 
are usually taken. Specifically, in estimating weak lensing 
masses one often bins galaxies in annuli, and estimates the 
mean shear in an annulus by dividing by the number of galax- 
ies found within the annulus, in other words, adopting a local 
normalization. Peak finders, on the other hand, often estimate 
the shear signal within a region by simply filtering the shear 
map which the galaxies sample, without necessarily dividing 
by the number of galaxies within the filtered region. In other 
words, they employ a global normalization, one that does not 
depend on the local galaxy density. In these papers, we ex- 
plore how each of these choices of normalization can impact 
the statistics and systematics of both peak finding and mass 
estimation. 

This paper is organized as follows. In section 2 we in- 
troduce the weak lensing estimator and calculate the statis- 
tical uncertainties of cluster mass obtained from this estima- 
tor. Various systematic errors are addressed in section 3. We 
study the impact of these mass estimates on the cosmological 
parameter constraints in section 4. Finally, we conclude in 
section 5. 

We adopt a fiducial flat ACDM cosmology with Q m = 0.28, 
Q A = 0.72, h = 0.7, AC = 4.54 x 10" 5 , and n s = 0.96. All 
masses are M2oo»S and distances are physical distances in ei- 
ther kpc or Mpc (as opposed to h~ l kpc or /r'Mpc). Finally, 
we set the lensing bias parameter q = 1.5 (see below for de- 
tails). 

2. MASS CALIBRATION IN LARGE OPTICAL SURVEYS VIA 
CLUSTER STACKING 

2.1. The Weak Lensing Shear Estimator 
8 i.e. 200 overdensity with respect to the mean matter density. 



We consider a shear weak lensing estimator of the form 

i 

Here, the sum is over all galaxies within some annulus of ra- 
dius R, e, is the ellipticity of galaxy i, h is the mean density of 
galaxies, and A is the area of the annulus. Expert readers will 
immediately discern that the estimator AE is different from 
the standard weak lensing shear estimator: 

i 

where N is the total number of source galaxies within the an- 
nulus. 

For the next few sections, we will focus primarily on AE 

rather than AS . This is primarily for convenience: the fact 

that the estimator AE takes the form AE = xjy where x and 
y are correlated implies that computing its mean and variance 
requires some additional algebraic gymnastics that we need 

not worry about when considering AE. Consequently, we 

have opted to illustrate our discussion with AE first, and then 

treat the more complicated case of AE . SectionOdiscusses 
how the two estimator s co mpare in terms of statistical preci- 
sion, while in section [3~4l we demonstrate that the two esti- 
mators have significantly different systematics with respect to 
foreground contamination of the source galaxy population. 

We wish to estimate the mean and variance of AE. In the 
interest of simplicity, we assume all sources are at the same 
redshift z and set E c (z,) = E c . It is also useful to rewrite our 
estimator as follows: first, we pixelize the sky behind the lens 
in pixels of area Ail, and define N, as the number of galaxies 
in pixel i. One might naively expect 

Ni = nMl{l + 5i), (3) 

where the (5, is the overdensity of matter in pixel i. In practice, 
however, lensing modifies the o bserved source galaxy density 
such that (ISchmidt et al.ll2009bh 

N i = n&VLtf l2 (l+8,) (4) 

where /i, is the magnification evaluated at pixel i, and q is 
a number that characterizes how the source density changes 
due to gra vitational lensing. For a DES-like survey, we expect 
q 1 -2 ( ISchmidt et al.ll2009bl) . Defining the filter function 
Wi such that = 1 only when a pixel is within the annulus of 
interest, and assuming the pixels are small enough that Ni is 

always either zero or one, we rewrite AE as 

AS = |^AO /2 ( 1 + SiWtet (5) 

i 

where the sum is now over all pixels. This is the equation we 
will use to derive the expectation value and variance of AE. 

2.2. Mean and Variance of AE 

The expectation value of AE follows from direct substitu- 
tion of the expectation values of e, and 5,. Specifically, 

(6) 

(Si}=0 (7) 
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where g is the reduced shear, g = 7/(1 — k). We assume now 
that we work in the thin annulus approximation so that the 
variation of fi and g within an annulus is negligible (though 
see also section 13731 . and find 



AS) = S c /^ 2 g«S f7 



(8) 



where all lensing quantities are to be evaluated at the annulus 
radius of interest. The approximate equality is obtained by 
expanding to leading order in the lensing quantities. 

To compute the variance, we assume source galaxies are 
unclustered, so that N, is Poisson. We have then 

Cov(N h Nj) = 5 y Var(iVj) = S u (Nt) = 6 tj fi 9/z nACl, (9) 



and therefore 



Cov(N h Nj) 

(Nt)(Nj) 



'V /2 «An' 



(10) 



In appendix [A] we consider the additional contribution to 
the noise due to clustering of the source population and we 
demonstrate that it is sub-dominant at high masses, and neg- 
ligible at low masses, justifying the assumption above. We 
further assume the variance in the ellipticity of galaxies is 
dominated by shape noise, 



{eiej)=gigj + Sij 



(11) 



the factor of two coming from the fact that there are two inde- 
pendent ellipticity components and we use only one of them. 
For quantitative purposes, we adopt a e = 0.3 as our fiducial 
value for the amplitude of the shape noise term. Direct sub- 
stitution results in 



Var(AS) = S 2 ^— 
nA 



2 1 7 



2hA' (12) 



where again the approximate equality gives the leading order 
term in a series expansion of the lensing quantities. 

The signal-to-noise of our estimator AS is 



m) l/2 i 



(nA) 1 / 2 



1 



(g 2 + ( r 2 /2) 1/2 (o2/2)V2 



(13) 



In a DES-like survey, we expect q w 1 -2, in which case the 
lensing bias correction \i q l A results in a modest increase of the 
signal-to-noise of our measurement. 

2.3. Survey Assumptions, Fiducial Model, and the Fisher 

Matrix 

We wish to estimate the precision with which stacked shear 
weak lensing experiments can constrain halo mass and con- 
centration in a DES-like survey. We adopt a survey area Q, = 
5,000 deg 2 , a fiducial galaxy density h= lOgalaxies/arcmin 2 , 
and a source redshift distribution 



f(z)<xz m exp[-(z/z*f] 



(14) 



with z* = 0.5, m = 2, and /3 = 1.4, as appropriate for a DES- 
like survey. The fraction of galaxies above a given redshift z 
is given by 



Fbg(z) = 



Io°°dz'f(z'y 



(15) 



Given a lens redshift zl, the (approximate) effective source 
density is simply fiFbgizi)- For a characteristic lens red- 
shift zl ~ 0.5, the corresponding source density is « 
8 galaxies/deg 2 . 

We further assume clusters are binned in narrow redshift 
slices z = z ± 0.05, and are logarithmically binned in mass in 
bins of width ±Alog 10 M = 0.1, corresponding to 5 bins per 
decade in mass. The effective source density for a given mass 
and redshift bin is therefore 



n eS = 2nF bg (z) 



dn 



A In M AV 



(16) 



where the factor of 2 arises from the total width of the mass 
bin 2 A InM, AV is the survey volume enclosed by the redshift 
slice z ± Az, 



AV = (1 +z) 2 D\(z)ncH-\z)2Az, 



(17) 



Da being the angular diameter distance, and Fb g (z) is the frac- 
tion of source galaxies with redshift higher than z- When es- 
timating all lensing properties, we will further assume that all 
source galaxies behind the clusters reside at a single source 
redshift z s equal to the mean source redshift of the sources be- 
hind th e cluster. The mass function dn jd InM is computed us- 
ing the lTinker et alJ (120081) mass function in our fiducial cos- 
mology. As mentioned in the introduction, for the purposes 
of com puting the impact of lensing bias we always assume 
2=1.5 (lSchmidtetalJl2009bl) . 

Note that we have not taken into account the effects of finite 
mass bins here: since the scale radius of halos evolves with 
mass, the stacked shear profile within [InM- A InM; InM + 
A InM] is not equal to the profile of an NFW halo with mass 
M. While this needs to be taken into account when fitting 
actual data, it is not of direct relevance to our Fisher fore- 
cast, so we neglect this effect here. Furthermore, it should 
be noted that, observationally, cluster stacks are made by bin- 
ning in an observable X, whereas the problem we have laid 
out here assumes that the clusters are binned in mass. While 
this obviously affects the distribution of cluster masses within 
a stack, we do not expect the precision of the corresponding 
measurement to be particularly sensitive to said distribution. 
Indeed, roughly speaking the total uncertainty in the mass per 
cluster is just the sum in quadrature of the measurement er- 
ror with the intrinsic scatter. For a typical cluster, the latter is 
smaller than the former, especially at low masses where most 
of our cluster sample resides, and therefore our approximation 
is valid. Given our expectations and the fact that our Fisher 
Matrix forecast ought to be interpreted as a rough estimate of 
the precision of these type of analysis, we have ignored this 
(mass proxy-dependent) complication. 

Assuming sources are uncorrected and that the redshift 
slices are narrow enough that halos are non-overlapping, the 

estimators AS for different radial bins are uncorrelated. In 
this limit, the Fisher matrix for a weak lensing shear experi- 
ment simplifies tc0 



F a b 



l d(AS a )d(AS Q 



(18) 



Var(AS a ) dp a dp b 

9 In the interest of simplicity, in equation [T8l we have neglected the mod- 
icum amount of information in the small dependence of Var(AS) on the 
model parameters p. We do not expect our results to be sensitive to this 
detail. Moreover, in practice we expect the covariance matrix for the mea- 
surements will be estimated using a jackknife method, which will erase any 
information in the variation of the covariance matrix with p. 
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where the sum is over all radial bins a and p is the vector of 
parameters of interest. Throughout, we assume logarithmic 
radial binning with bins of width ±Alog 10 /? = 0.02. More- 
over, we will only add bins over the radial range 0.1 Mpc to 
2 Mpc. We will further impose the condition that the clus- 
ter magnification must be smaller than five (i.e. < 5) for 
a radial bin to be included in our computation. This ensures 
that for those systems where the Einstein radius extends past 
the 0.1 Mpc minimum radius, we do not include information 
from sources that are strongly lensed. Our results are robust 
to making our bins narrower, but do depend in de tail o n the 
radial range assumed. For a discussion, see section [372] 

We model the mass distributio n of cluster stacks as a 
Navarro, Frenk, and White (NFW, Navarro et ail 1 19961) pro- 
file. For simplicity, we have held the concentration parameter 
of our clusters fixed to c = 5. Changes in concentration within 
a stack do not qualitatively change any of our results, and have 
only a modest impact on the quantitative results. All of our 
conclusions are robust to the choice of concentration param- 
eter. We compute the c onvergence and shear of N FW pro- 
files using the form ulae in lWright & Brainer dl(l2000T) ( see also 
lBartelmannlll996L and appendix of iSchmidt & Rozol ( 120 101) ). 

2.4. Results 

Figure Q] shows the predicted statistical uncertainty in the 
weak lensing mass and concentration of our cluster stacks as a 
function of the mean cluster mass of the stack. Different lines 
correspond to different lens redshift. From top to bottom, the 
solid lines assume lens redshifts ziens = 0.1, 0.2, and 0.4. From 
bottom to top, the dashed lines assume z\ em = 0.6, 0.9, and 
1.2. The vertical dotted line is a rough estimate of the ex- 
pected mass threshold for optical selection. Qualitatively, the 
precision of weak lensing measurements increases from z = 
to z ~ 0.2, reflecting the increasing number of lenses due to 
increased survey volume, as well as the improvement on the 
lensing efficiency of the lenses. Between z ~ 0.2 and z ~ 0.6, 
the precision of the weak lensing measurements is roughly 
redshift independent: even though survey volume continues 
to increase, it is now offset by a diminishing source density. 
For redshifts z > 0.6, both the lensing efficiency and the effec- 
tive source density decrease quickly with increasing lens red- 
shift, so weak lensing masses begin to worsen. We also note 
that at all redshifts there is a mass scale at which the errors 
blow up, reflecting the exponential drop off in the halo mass 
function. The typical precision with which the mean mass 
of galaxy clusters may be estimated at moderate redshifts is 
roughly 2%. 

2.5. On the Choice of Estimator 



In section ETI we defined the estimator AS differently from 
the more commonly used weak lensing shear estimator: 



(19) 



(20) 



Pixelizing the sky, the estimator AS can be written as 

^/ = - E,.Ai^f /2 (i+&)w;- g< - 

where the sum is now over all pixels. 

A difficulty with AS is now immediately apparent. We see 
that AS takes the form AS = x/y, and in general one has 
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FIG. 1. — The predicted 1<t statistical error in the log-mass (top) and the 
concentration (bottom) of galaxy clusters in a DES-like survey as measured 
through cluster stacking. The solid lines illustrate the trends with redshift for 
Z < 0.4, while the dashed lines illustrate the trends with redshift for z > 0.4. 
From top to bottom (in either plot) the lens redshifts of the solid curves are 
Hens =0.1, 0.2, 0.4, whereas the redshift of the dashed curves is, from bottom 
to top, ziens =0.6, 0.9, and 1.2. The typical precision with which the mean 
mass of galaxy clusters may be estimated at moderate redshifts is roughly 
2%. 



AS J = (x/y) =/ (x) I (y). One way to address this problem 

is to write y=(y) + Sy, and expand the denominator in a power 
series using the binomial expansion. This results in a power 
series expansion of all quantities in terms of 1 /(nA) (see Ap- 
pendix [Bj. Doing so up to second order in 8 to compute the 

mean of AS we find 



AS > = S c 



(21) 



The quantities in the angular brackets in the right hand side of 
the equation are to be averaged over the radial bin of interest. 
Adopting the thin annulus approximation, we can pull out the 

//f term out of the sums in equation|20]to obtain 

£,.Aft(l + <5,> ; 



AS = Y, C '- 



(22) 



In this limit, the magnification term drops entirely out of the 
equations, and one can easily solve for the mean and variance 

of AS using the power series approach advocated earlier. We 



On Stacked Weak Lensing 



5 



find (see AppendixIBl 



AS ) = E c g 



Var(AS ) = £? 



1 



9/ 2 «A 2 



(23) 



(24) 



Note the magnification still appears in the variance since the 
effective source density is [i q l 2 n. 

Using the above equations, we repeat the Fisher matrix 
experiment that we carried out for AS. We find that these 
two estimators recover the mean cluster mass of a stack with 
nearly identical precision, with \om-om 1 | 0.2% where <tm 
and <tm' are the forecasted uncertainties in (InM) for a stack 

using the estimators AS and AS respectively. 

Before we end, we would like to put a word of warning con- 
cerning the validity of our results for the estimator AS . Our 
derivation in this section deals with the covariance between 
the numerator and denominator of AS through a power series 
expansion, where the expansion parameter is 1 / nA, the ex- 
pected number of galaxies within the annulus of interest. Con- 
sequently, for our analysis to be applicable one must choose 
radial bins that are large enough for ilA ^> 1. On the other 
hand, we are also using the thin annulus approximation, so A 
cannot be arbitrarily large. While we expect to be able to si- 
multaneously satisfy both of these constraints in stacked weak 
lensing analysis, this is generally not possible for individual 
clusters. For instance, in our fiducial model, a relatively broad 
radial bin R e [0.2 Mpc,0.4 Mpc] would only contain w 25 
galaxies, so the higher-order terms in the 1 / ilA power series 
expansion can in principle introduce w 5% level corrections 
to the signal in the inner most radial bins. Here, we do not 
address this additional difficulty since our primary interest is 
stacked weak lensing. The finite width bin corrections that 
we treat below concern corrections to the leading order term 
which occur even when nA 3> 1 . 

3. SYSTEMATICS 
3.1. Lensing Bias Corrections 
By lensing bias corrections we refer to the terms that scale 
as [i q l N for some value of N in both the mean a nd v ariance 
of our shear estimators. As discussed in section 12.51 in the 



thin annulus approximation the mean of AS is independent 
of these lensing bias terms, so we expect lensing bias correc- 
tions to be smallH The same is not true of AS, and therefore 
ignoring these corrections can potentially lead to systematic 
biases. That said, it is worth emphasizing that incorporating 
these corrections into the analysis is not difficult, so having 
to include such corrections is not any sort of limiting system- 
atic. Our goal here is simply to determine whether doing so is 
a necessary step in the analysis of future data sets. 

To address this question, we use the Fi s her m atrix for- 
malism set out in the ap pendix of IWu et ail (120081) (see also 
iHuterer & Linderl 120071) to estimate the mass and concen- 
tration parameters that would be recovered from the data if 

10 In principle, the fact that lensing bias does impact the covariance matrix 
of AS could lead to biasing in halo mass and concentration. In practice, 

however, the error bars employed when relying on AS are based on the 
Poisson variance of the obsen'ed number density, which includes the appro- 
priate lensing bias. Thus, one does not expect AS to be affected by lensing 
bias. 
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FIG. 2. — 68% likelihood contours obtained when measuring the mass and 
concentration of a halo stack of mass M = 5 X 10 14 Mq at redshift z = 0.3 us- 
ing a DES-like survey. The solid contour includes the corrections for lensing 
bias, while the dashed line ignores these corrections. The best fit value for 
the former case coincides with the input parameters, marked by a diamond, 
while ignoring these corrections leads to biased expectation values, shown 
here with a triangle. This behavior is generic (see Figuref5]for details). 

lensing bias corrections are ignored, and then compare those 
results to those derived when properly including the lens- 
ing bias terms in the analysis. The basic idea is illustrated 
in Figure [2] The solid and dashed ellipses represent the 
68% confidence contours when including (solid) or ignoring 
(dashed) lensing bias corrections for a halo stack at z = 0.3 
and M = 5 x 10 14 Mq. The true halo mass and concentration 
is marked by the diamond, while the triangle marks the best 
fit halo mass and concentration obtained when ignoring lens- 
ing bias corrections. We can see that the best fit values when 
ignoring lensing bias effects are well beyond the 68% error 
ellipse of the experiment, demonstrating that lensing bias is a 

significant correction when using the AS estimator. 

The extent to which lensing bias is important relative to sta- 
tistical uncertainties depends on both halo mass and redshift. 
In particular, lensing bias is more important at high masses, 
and at redshifts z ~ 0.2-0.5, for which lensing efficiency is 
high and statistical errors are small due to high source den- 
sities. This is illustrated in Figure [3] where we show the 
ratio between the systematic bias in mass and concentration 
incurred by ignoring lensing bias effects, relative to the cor- 
responding statistical uncertainty. We see that our projected 
bias in the halo mass can be twice as large as the correspond- 
ing statistical uncertainty, while the corresponding errors in 
concentration can be even larger than 10 times the statistical 
error. Thus, inclusion of lensing bias and reduced shear cor- 
rections are necessary for weak lensing mass calibrations in a 
DES-like survey. 

One aspect of our results that is very surprising is that lens- 
ing bias can lead to systematic underestimates of the mass, 
despite the fact that lensing bias leads to an apparent in- 
crease of the lensing signal. The reason we find that masses 
are underestimated when ignoring lensing bias is due to the 
anti-correlation between halo mass and concentration. Since 
lensing bias dramatically boosts the recovered concentration, 
the recovered mass must necessarily go down to compen- 
sate. This statement, however, is very much dependent on 
the range of radii considered in the calculation. In particular, 
as one moves out the maximum radius employed in the anal- 
ysis the anti-correlation between mass and concentration de- 
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FIG. 3. — Ratio of the systematic bias in the log-mass AlnM or concentra- 
tion Ac incurred when ignoring lensing bias and reduced shear errors, relative 
to the statistical uncertainty <T[ nM or a c in these quantities (see FigurefT). 



creases, and lensing bias beg ins t o produce positively biased 
halo masses (see also section [3T2T >. 

Finally, our r e sults may appear to contradict those in 
ISchmidt & Rozol (120101) . who find lensing bias results in an 
overestimate of weak lensing masses despite the fact that they 
estimate masses using filters of comparable radius to those 
employed here. The reason for this apparent discrepancy is 
that in this analysis we ha ve fit for both halo mass and con- 
centration, whereas in the ISchmidt & Rozol (120101) analysis, 
the concentration parameter is held fixed at its fiducial value. 
In that case, the apparent increase in the matter density in the 
inner regions of a cluster due to lensing bias necessarily goes 
into halo mass rather than concentration, leading to a positive 
mass bias. In practice, because cluster stacks have abundant 
signal, one fits for both mass and concentration, so the obser- 
vationally relevant case for cluster stacks is that considered 
here. For individual weak lensing mass estimates in the low 
signal-to-noise regime, the results of I Schmidt & Rozol (120101) 
are the relevant ones. 

3.2. The Relative Information Content of Different Scales 

In light of the results from the previous section it is worth 
investigating the contribution to the Fisher matrix from each 
individual radial bin as a function of radius. The information 
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FIG. 4. — Mass (solid) and concentration (dashed) sensitivity as a function 
of radius for a 10 15 Mq halo at redshift z = 0.6. The sensitivity functions 
characterize the relative information content of various scales. Dotted curves 
are the two sensitivity functions we obtain when neglecting lensing bias cor- 
rections. Not surprisingly, the scales that dominate the concentration and 
mass information content are a few hundred kpc and a few Mpc respectively. 
In both cases, lensing bias shifts the sensitivity functions towards smaller 
scales, as we would expect. 



contributed by radial bin R a is 



Fab(Ra) '■ 



1 



Var(AS Q ) dp a 



dpb 



(25) 



Of course, this quantity explicitly depends on the width of the 
radial bins employed. We therefore define the relative sensi- 
tivity s a b via 

s ah (R a )= Fa f fi, x (26) 

assuming logarithmic radial bins. Note that, by construction, 
Sab is ver y nearly independent of the assumed width of the 
radial bins and the assumed source density. Of particular in- 
terest to us are the mass sensitivity s mm and concentration sen- 
sitivity s cc , which characterize the relative mass and concen- 
tration information content of various scales. 

Figure 0] illustrates the mass (solid) and concentration 
(dashed) sensitivity functions for a 10 15 Mq halo at redshift 
Z = 0.6. Not surprisingly, the mass information content is 
dominated by scales of order a few Mpc, while the concen- 
tration information content is dominated by scales of order 
a few hundred kpc. The large overlap of the two sensitiv- 
ity functions results in the anti-correlation between halo mass 
and concentration noted earlier. The dotted lines in the fig- 
ure correspond to the sensitivity functions in the absence of 
lensing bias effects (i.e. setting q = 0). We see that lensing 
bias shifts the sensitivity functions to smaller scales, effec- 
tively up-weighting them, as we would expect. The corre- 
sponding sensitivity functions for the estimator AE are es- 
sentially identical to those of AE in the absence of lensing 
bias corrections. 

In light of Figure [4] it is easy to understand why our re- 
sults were sensitive to the adopted radial range. It is evident 
from the figure that both the minimum and maximum radius 
will impact our results. The minimum radius chosen will di- 
rectly affect the precision of the concentration measurement, 
which in turn propagates to the cluster mass via the mass- 
concentration correlation. Likewise, it is easy to see why our 



On Stacked Weak Lensing 



7 




0.00 0.02 0.04 0.06 0.08 0.10 
Logarithmic Width of Radiol Bin (±A log 10 R) 

FIG. 5. — Bias in the recovered weak lensing masses using the thin-annulus 

approximation for predicting the expectation value of AE as a function of 

the width of the annulus used to estimate the shear profile. We assume AS 
is estimated in annuli of width A defined via [\nR c - A, lni? c + A]. The solid 
line assumes the thin annulus approximation at R = R c , while the dashed line 
assumes the thin annulus approximation at the mean radius of the annulus 
(R). This bias is nearly independent of halo mass, concentration, redshift, 
and source density. 



Our results are shown in Figure [5] The largest bin width 
we consider is roughly five bins per decade in radius, corre- 
sponding to ±Alog 1() /? = 0.1. For such a bin width, there is 
a negative mass bias of « 5% assuming the thin annulus ap- 
proximation at R = R c (solid line), and a rj 2% positive bias 
assuming R = (R) (dashed line). This bias is nearly indepen- 
dent of halo mass, concentration, redshift, or source galaxy 
density. Thus, it is important to account for finite bin width 
effects when estimating cluster masses, particularly for the 
purposes of calibrating cluster scaling relations PI 

The origin of the the above bias is simple to understand: 
consider first the case R = R c . It is evident that in this case our 
annulus contains more area beyond R = R C than below R = R c . 
Consequently, we expect a negative bias, since the outer re- 
gion of the annulus is being up-weighted relative to the inner 
region. Setting R = (R) evidently overcompensates for this ef- 
fect. While one can in principle imagine finding the specific 
value of R for which there is no bias, this value will depend 
in detail on cosmology, source galaxy density (due to mag- 
nification bias corrections), etc., so it is simplest to just take 
into account the finite bin-width correction exactly. Alterna- 
tively, when performing stacked weak lensing measurements 
one may adopt very narrow radial bins, in which case these 
corrections become negligible. 



results are sensitive to the maximum radius employed: scales 
as large as tens of Mpc still contain non-negligible amounts of 
information pertaining to the halo mass. Of course, in practice 
such sc ales show devia tions from a simple projected NFW 
profile (iHoekstral (120031) ). both because of possible deviations 
from NFW in the 1-halo term, and because of the appearance 
of a 2-halo term. We defer an investigation of this problem to 
future work. Here, we simply note that we expect that there is 
significant additional information in the weak lensing signal 
up to scales as large as several tens of Mpc. 

3.3. The Thin Annulus Approximation 

When performing our forecasts for both AE and AS we 
adopted the thin annulus approximation, in which one ignores 
variations of the convergence, shear, and magnification fields 
within an annulus. When using finite radial bin widths, how- 
ever, there is a slight ambiguity on what radius should one em- 
ploy when performing the thin annulus approximation. Two 
obvious choices are the central radius R c where the radial bins 
are defined by \nR c ± A In/?, and the mean radius of the an- 
nulus 



(R): 



2Rl 



-Rl 



3R 2 



Rl 



(27) 



The question then becomes, how thin must an annulus be for 
the thin annulus approximation to hold for either of these two 
choices of radii? 

To address this question, we have again relied on a Fisher 
matrix approach to estimate the systematic error in the log- 
mass of a galaxy cluster when the mass is estimated using the 
thin-annulus approximation, as a function of the width of the 

annulus used to estimate AE or AE . For specificity, we will 

focus here on AE , but very similar conclusions hold for AE. 
We emphasize that, just as with lensing bias corrections, one 
can in principle just include these finite bin width effects in the 
fits. The analysis described here simply addresses whether it 
is necessary to do so. 



3.4. Foreground and Cluster Member Dilution 

One important systematic when estimating the mass of 
galaxy clusters throu gh weak lensing is the effect known 
as me mber dilution (iBernardeaul ( 11998 ): Medezinski et al. 
( 120071) ). That is, if cluster members or other foreground 
galaxies are included in the source population, they will re- 
duce the observed shear and cause us to underestimate the 
cluster mass. Let us treat this problem explicitly: the total 
(angular) galaxy density around a cluster at redshift z can be 
written as a sum 



rc = nbg + «fg + "ci 



(28) 



where «bg is the background galaxy density field, «f g is the 
foreground galaxy density, and n c \ is the density of cluster 
galaxies. When estimating weak lensing masses, one adopts 
photometric redshift cuts to preferentially select background 
galaxies. Let r^ % , rf g , and r c \ be the probabilities of passing 
the photometric redshift cuts for a galaxy in the background, 
foreground, and at the lens redshift, respectively. A perfect 
photometric redshift selection would result in rf g = 0, r c \ = 0, 
and Tb g = 1, so the value of these parameters characterizes the 
actual redshift distribution of the source sample in the pres- 
ence of imperfect photometric redshift selection. We wish 
to consider the impact of uncertainties in these parameters on 

AE and AE . For simplicity, we will work in the thin annulus 
limit throughout. 

We begin by considering AE, which we now write as 



AE = ■ 



nAF 



£AQ[ 



r bg n bg e,. fc + r fg 77f g e,. £ 



nAF ^ 



(29) 



11 The small "steps" that appear in Figure|5]are real. To estimate cluster 
masses, we use the scales 0.1 Mpc < R < 2 Mpc. This means that when we 
vary the bin width, there will be discrete jumps in the number of bins, which 
leads to the "stepping" observed in the plot. 
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FIG. 6. — Top Panel: The ratio n c i/n of the projected density of cluster 
galaxies to the mean galaxy density of the survey for a halo of mass M = 
10 15 Mq for a variety of redshifts. From bottom to top, the redshift of the 
solid curves are z =0.1, 0.15, 0.3, and 0.6. From top to bottom, the redshift of 
the dashed lines are z =0.9, 1.2, and 1.5. The redshift slice 0.2 < z < 1.0, the 
ratio n c i jn can be larger than unity in the inner few hundred kpc for massive 
halos. 

Here, F is the fraction of galaxies with redshift higher than the 
redshift of the cluster (according to the observed redshift dis- 
tribution), while n is the mean galaxy density over the whole 
survey. Hence, averaging over the whole survey, 



We will again ignore correction terms to the mean that go as 




(n hg )=Fn; (n fg ) = (1 -F)n. 



(30) 



Further, e fg , e bg , and e cl represent the foreground, background, 
and cluster ellipticity fields. The first two moments of e bg sat- 
isfy equations l6l and [TT1 while for the foreground and cluster 
ellipticity fields we set 

(31) 
(32) 



(*,■)= 
{e i e j ) = -o 2 e 8 i j. 



In practice, cluster member galaxies are known to be radially 
aligned (i.e., (ef 1 ) =^0), but the extent of the radial ellipticities 
is less well known. We defer this additional source of error to 
future work. Upon plugging in, we find 



AS ) = r bg E, 



(33) 



implying that AE is sensitive to only one of the three system- 
atic parameters considered here. The relative systematic error 
due to an uncertainty <5r bg in rb g is simply 



cr sys (ln AS) = — - = Sr bg 



(34) 



assuming our fiducial value rb g = 1. 
Let us now turn our attention to AE , which takes the form 

A" 



where 



AE =E C - 

y 



rb s n bg e^ s + r {g n {g e? + r c \n c \ef 



i 

y = y^AO [r bg n bg + r fg n fg + r ci n d ] 



(35) 

(36) 
(37) 



1 /(«A), and set ( AE 



AE ) = E e g 



1 + 



= (x) I (y). We then arrive at 
r fg l-F r d (« d ) 1 



(38) 



r bg )Jfil 2 F rb g n ^fll 2 F 
Here, (n c \) is the mean number density of cluster member 
galaxies in the stack considered. We find that ^ AE \ is sensi- 
tive to all three of our systematics parameters rf g , r bg , and r c \. 
While these parameters only appear as two independent ratios 
— rf g /rb g and r c i/rb g — in practice we still need to estimate 
all three parameters independently. 

But just how sensitive is AE to systematic uncertainties? 

As for AE, we compute the relative systematic error due to a 
small uncertainty in our systematics parameters. Linearizing, 
we find 



CT sys (lnAE') = -<5r ft 



l-F 



■Sr c i 



(«d> 1 



(39) 



yfl^F h ^l 2 F' 

Comparing the above expression to equation[3H it is clear that 

AE is more robust to systematics if the coefficient of Srf g and 
Sr c \ are smaller than unity. Inspecting the above expression 
we notice two important things: first, the coefficients scale 

as l/F, so at large redshifts, AE is guaranteed to become 

more robust than AE . Second, the coefficient for Sr c \ scales 
as n c \/n (dilution effect), which we expect to be larger than 
unity in the cores of clusters, but insignificant in the outskirts, 

so it is possible for AE to be more robust than AE in the 
cores of clusters, but the opposite be true in the outskirts. 

Figure [6] shows the ratio n c i/« for a halo of mass 10 15 M© 
as a function of radius for a variety of redshifts. The ratio has 
been estimated using the model detailed in Appendix ICl As 
we expected, this ratio is larger than unity in the inner few 
hundred kpc for 0.2 < z < 1 .2. Note, however, that this is true 
for very massive halos. If we assume instead M = 10 14 , this 
ratio barely reaches unity at R = 100 kpc. Note that the simple 
model for cluster members used here should only be taken as 
an order-of-magnitude estimate. 

Figure [7] shows the relative systematic error for both the 

AE and AE estimators assuming a systematic uncertainty 
S rbg = 5 rfg = <5r c i =1% (one can easily scale the values in Fig- 
ure |7]to other valus of Sr). For the estimator AE, the corre- 
sponding systematic error is constant, and mass and redshift 

independent. The same is not true of AE . We find that for 

massive halos, the error in AE is typically dominated by di- 
lution of the lensing signal by cluster member galaxies. This 

error becomes larger than the expected statistical error in AE 
at z ~ 0.3, though it affects only the innermost few hundred 
kpc scales. Moreover, this error is explicitly mass dependent, 
with less massive halos benefiting from smaller systematic 
uncertainties. In all cases, however, for z > 0.8, the l/F de- 
pendence of the systematic error on AE makes this estimator 

less robust than AE. 

From the above discussion alone, it is unclear if either of 
the two estimators we have considered here is superior to 
the other. At very high masses and high redshift, AE is 
likely to be preferable to AE , but this may be reversed as 
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FIG. 7. — Relative systematic error for both AS (dotted) and AS (solid 
and dashed) as a function of radius for a halo of mass M= 10 15 Mg. The solid 

curves show the error on AS due to a 1% systematic uncertainty Srf„ on the 
fraction of foreground galaxies included in the source galaxy population (the 
curves for z = 0.1,0.3 are below the minimum error plotted). Dashed curves 

are the relative systematic error on AS due to a 1% error on r c i, the fraction 
of cluster galaxies mistakenly included in the source population. Finally, the 
dotted curves are the (redshift-independent) relative error on AS due to a 1% 
systematic error 8r\,„ in the estimated fraction of background galaxies in the 
source population. 

one moves towards lower redshifts and/or masses. Further- 
more, choosing between AS and AS also amounts to choos- 
ing between a spatially constant, mass-independent but likely 

redshift-dependent systematic in the case of AS, and a sig- 
nificantly mass- and radius-dependent systematic in case of 

AS . The cluster member dilution affecting AS is frequently 
corrected for by multiplying the shear profile measured with 

AS by the observed source galaxy-cluster correlatio n func- 
tion (e.g., "boost factor" in iMandelbaum et alJ 120051) . This 

is essentially equivalent to using the estimator AS from the 
start. Finally, the answer as to which estimator is preferable 
will also depend on which of the systematic uncertainties dis- 
cussed here — r\, g , rf g , and r c \ — can be best controlled. 

3.5. Other Systematics 

In addition to the systematics discussed above, there are 
additional sources of systematic uncertainty that can impact 
our results. A few examples of such sources of systematic 
uncertainty are: 

• Shear Systematics: if shear is mis-estimated, this sys- 
tematic will of course be carried over to the estimated 
cluster mass. 

• Miscentering: if clusters are miscentered, this can have 
a dramatic impact on the expectation va lue of the weak 
lensing signal in the cores of clusters dJohnston et al.l 
l2007bl) . This systematic will almost certainly dominate 
the uncertainty with which the concentration of galaxy 
clusters can be measured within cluster stacks, but its 
impact on cluster masse s can be significantly red uced 
through careful analysis (Mandelb aum et alj|2010l) . 

• Source obscuration by foreground galaxies: occasion- 
ally, background galaxies will be perfectly aligned with 
foreground galaxies, and therefore the latter cannot be 



included in the weak lensing shear signal estimation. 
This reduces the effective area of the annuli used to esti- 
mate shear, which in turn impacts the expected number 
of galaxies in the annulus, thereby impacting globally 
normalized estimators. 

• Photometric redshift errors: the strength of the shear 
signal depends on the redshift of the source galaxy 
under consideration, which must be estimated based 
on photometric data. Consequently, scatter and catas- 
trophic errors in photometric redshift estimates may 
have a significant i mpact on shear mass cali bration ex- 
periments (see e.g. [Mandelba um et al.ll2008b . 

In light of this discussion, we reiterate that determining 

whether AS is superior to AS or vice-versa will require 
empirical investigation, and that due to these systematics, the 
forecasted precison for a DES-like survey is best interpreted 
as the lower-limit for what will actually be realized. That said, 
until one estimator may be conclusively shown to be superior 
to the other, estimating cluster masses using both estimators 
should allow one to estimate the level of systematic uncer- 
tainty introduced by source galaxy selection, which should be 
of tremendous utility in upcoming photometric surveys. 

4. THE IMPACT OF STACKED WEAK LENSING ON CLUSTER 
ABUNDANCE EXPERIMENTS 

4. 1. Fisher Matrix and The Fiducial Model 

We now consider whether the statistical precision of the 
mean weak lensing masses recovered from stacked weak lens- 
ing is sufficient to significantly improve cosmological con- 
straints in a DES-like survey relative to the self-calibration 
expectation with Planck priors. We address this question by 
once again relying on the Fisher matrix formalism. Specifi- 
cally, we set the Fisher matrix for our experiment to the sum 
of the standard self-calibration result plus an additional con- 
tribution due to the stacked weak lensing data, 



Ftotal - Fpi anc k + Fsc + F\ 



WL- 



(40) 



The Planck priors Fisher matrix is that provided by Hu and 
Ma (private communication). For the self-calibrat ion fisher 
matrix , we use the formalis m described in detail in IWu et all 
( 120081) and lWu etail (120101) . Briefly, we assume that the sur- 
vey area 5,000 deg 2 is divided into 500 patches of 10 deg 2 
each. The observables in our experiment are the cluster counts 
in each of these 500 patches for each mass and redshift bin. 
We adopt the same binning as for the weak lensing analy- 
sis (five bins per decade in mass, and redshift slices of width 
±Az = 0.05). 

The model parameters we consider can be split into two 
categories, cosmological parameters, and the nuisance param- 
eters describing the observable-mass relation. The cosmo- 
logical parameters and their fiducial values are given in the 
introduction, except that we add as free parameters wq and 
w a . These describe the equation of state of the dark energy 
w = wo + w n ( 1 - a) as a function of the cosmic expansion fac- 
tor a. We wi ll be primarily interes ted in how the Dark Energy 
Task Force (lAlbrecht et al .1120061) figure of merit — defined 
as the the product of the eigenvalues of the w and w a Fisher 
matrix — changes upon inclusion of the weak lensing data. 

In order to describe the observable-mass relation, each 
cluster is assumed to have an observed "mass" M„/„ which 
represents its richness measurement. We assume the 
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observable-mass relation P(M obs \M true ,z) is log-normal, with 
the mean and variance of lnM b s assumed to scale linearly 
with InM and ln(l +z). We write 



(lnM ofc! )=lnMo + Q:Mln 



M 

M pivot 



+ a z \n(\+z) (41) 



M 



M pimt 



+ /3 z ln(l+z). (42) 



We set Mptvot = 7 x 10 13 MqQ The remaining set of parame- 
ters are allowed to vary, and all of our results are marginalized 
over these observable-mass parameters. We set the fiducial 
value of our free parameters to those of an unbiased estimator 
in InMobs, so that InMo = InM and % = and a z = 0. Further, 
for our fiducial model we assume no evolution of the scatter 
with mass or redshift, so that Pm = Pz = 0- We consider two 
possible values of the amplitude of the scatter <to = 0.2 and 
(To = 0.5, corresponding to 20% and 50% scatter in M obs at 
fixed M true . 

Our stacked weak lensing analysis allows us to introduce an 
additional set of observables in our analysis, namely the mean 
mass Mwl in bins of M obs . Since the mean mass estimates of 
the different bins are independent, the weak lensing Fisher 
matrix is simply 



Fwr. — 



E 

all bins 



1 d(lnM WL ) d(lnM WL ) 



'WL 



dpi 



d Pj 



(43) 



where a WL is the error in the mass estimated in section [2~41 To 
compute (In Mwl), we assume that weak lensing masses are 
unbiased, so that 



(M WL \M obs ,z) = (M true \M obs ,z) 



(44) 



If dn dV , . 

= T / dM trm dzM true — — — (4>) (45) 
N J dMtrue dz 

where <j> is the binning function in observed mass and redshift, 
(4>) is the effective binning function as a function of true halo 

mass, 



(<p\M mie ,z)= dM obs (j)(M abs , z)P{M obs \M tme , z) , (46) 



and N is the expected number of clusters, 

dn dV 



N= dM 1rue dz 



dMt rue dz 



(47) 



A more detailed discussion of how to add t his additional in - 
formation to cluster forecasts is presented in lWu et all (12010b . 

Before moving on, we point out that since the err or e sti- 
mates <jwl of the weak lensing masses from section 12.41 as- 
sumed clusters were binned according to their true masses, 
there is a small scatter-dependent correction to the predicted 
uncertainties. Given that we have n ot ta ken into account finite 
mass bin width effects (see section l2~3l l. we ignore these cor- 
rections in this section as well, and simply remind the reader 
that the corrections are expected to be small since shape noise 
is larger than the intrinsic scatter in the mass. 

4.2. Results 



Note InMo and lnM p ,,, or are degenerate, so we can fix one of them arbi- 
trarily without loss of generality. 



1 1 ~> 



n src - 2 n Fiduciol 

Fiducial 



n =0.5n 



src -u --"'Fidudol 



-3 -2-10 1 

log 10 (Variance Prior) 

FIG. 8. — The figure of merit for a DES-like cluster abundance experiment 
with stacked weak lensing mass calibration relative to its self-calibration ex- 
pectation, as function of the external prior Actq on the scatter in the mass- 
observable relation. Planck priors a re as sumed. The solid lines employ the 
forecasted errors arrived at in section fOl The dashed and dotted lines have a 
source density that is a factor of two lower (dashed) and higher (dotted) than 
that of our fiducial model. The long tick marks along the top axis correspond 
to a 5% prior on the scatter assuming cq = 0. 1 (solid), <ro = 0.2 (dashed), and 
Co = 0.5 (dotted). The relative improvement in the figure of merit is nearly 
independent of the adopted value of ctq, particularly for broad priors. 



Figure [8] shows the figure of merit of our DES-like survey 
including stacked weak lensing mass calibration relative to 
the figure of merit obtained without this additional source of 
data. We have allowed for the possibility of a prior on the 
scatter parameters with Actq = A0m = A/3 Z . The solid line 
assumes the weak lensing uncert ainty for each mass and red- 
shift bin <jwl estimated in section l2~4l whereas the dashed and 
dotted assumed a source density that is half and twice that of 
our fiducial model respectively. Assuming no scatter prior, 
we find that the improvement in the Dark Energy Task Force 
figure of merit is in the range 2.5 - 3.5. This is true for both 
the ao = 0.2 and the ao = 0.5 models: the relative improve- 
ment in the figure of merit is only weakly dependent on the 
assumed scatter, which is why we only show one set of curves 
in Figure[8](those for ao = 0.5). The fiducial value for the fig- 
ure of merit is, of course, different. We have FoMsc =19 for 
do = 0.2, while FoMsc = 15 for <to = 0.5. 

Additional priors on the scatter of the mass-observable re- 
lation may further increase the efficacy of stacked weak lens- 
ing mass calibration, leading to improvements in the figure of 
merit as large as a factor of 8. To do so, however the priors 
need to be very tight. Note that we have expressed these pri- 
ors in terms of the variance rather than the standard deviation. 
For reference, the long thin vertical tick marks a long the top 
x-axis in Figure [8] correspond to a 5% prior (i.e. Ado = 0.05 
on the scatter, assuming ao = 0.1 (solid), 00 = 0.2 (dashed), 
and (To = 0.5 (dotted). Note that while strictly speaking we 
are only showing the relative improvement for the ao = 0.5 
case, the curves for (t = 0.1 and ao = 0.2 closely track that for 
(To = 0.5. It is evident from the Figure that accurate knowledge 
of the scatter can lead to significant further improvement in 
the figure of merit of the experiment. 

We have also considered how biases in the weak lensing 
masses could impact cosmological parameter estimation in 
order to assess the level at which systematics need to be con- 
trolled. Not surprisingly, we find that « 2% biases in mass — 
which correspond tow la — result in « la biases in the in- 
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ferred cosmological parameters. Whether the recovered weak 
lensing masses can indeed be expected to be unbiased at the 
2% level in real data, however, remains to be seen. 

5. SUMMARY AND CONCLUSIONS 

Weak lensing shear profiles and weak lensing peak finding 
tend to utilize the shear signal in different ways, with shear 
profile mass calibration relying on locally normalized estima- 
tors (AS ), and peak finding relying on globally normalized 

estimators (AS). We have used both of these estimators to 
predict the precision to which the mean mass of galaxy clus- 
ter stacks in a DES-like survey can be measured. We find 
that for moderate redshift clusters (z < 0.6), the typical preci- 
sion achieved is « 2%, with the two types of estimators hav- 
ing nearly identical statistica l uncertainties. A companion pa- 
per, [Sch^dL&Roi3 (120101) . investigates similar issues in the 
case of weak lensing peak finding. 

We also considered three sources of systematic biases for 
these measurements. The first is lensing bias, which we find 

affects AS but not AS . For the former, we find that includ- 
ing lensing bias corrections to the expectation value of AS is 
necessary to avoid significant systematic biases in both halo 
mass and concentration. We emphasize however that lensing 
bias corrections can be easily incorporated into the data anal- 
ysis, so this is not a particularly worrisome systematic. Our 
results simply indicate that it is necessary to incorporate these 
corrections. 

The second source of systematic uncertainty we consid- 
ered are finite bin-width corrections, which affect both AS 

and AS . Again, these corrections can in principle be ex- 
plicitly included when analyzing data, and our investigation 
only addresses whether doing so is necessary for practical pur- 
poses. We find that for logarithmic bins ±Alog 10 /? < 0.04, 
the systematic bias in mass from finite bin-width corrections 
is less than 1%. If one uses 5 bins per decade in radius 
(±Alog 10 7? = 0.1), biases as large as 5% in mass are possi- 
ble. 

The final systematic we consider here is fluctuations in the 
number density of galaxies due to intrinsic clustering, which 



can affect the estimators through imperfect photometric red- 
shift selections. Remarkably, we find that AS and AS are 

affected by this systematic in very different ways: AS is af- 
fected by the fraction of background galaxies missed by pho- 
tometric redshift selection, whereas AS is affected by the 
fraction of foreground or cluster member galaxies that are in- 
cluded in the source population. Whether AS or AS is more 
robust depends on which of these systematics can be better 
controlled, a question which can only be empirically resolved. 
In either case, both estimators can be used to cross-check each 
other for the effects discussed here, due to their significantly 
different systematics. 

Finally, having estimated the precision with which the 
mean cluster mass of clusters stacks can be measured with 
a DES-like survey, we have investigated how this measure- 
ment would impact the cosmological parameter constraints of 
a DES-like cluster abundance experiment. We find that with 
our fiducial assumptions, the figure of merit of such an exper- 
iment improves by a factor of 2.5-3.5, with larger increases 
possible if priors on the scatter of the observable-mass re- 
lation can be derived from additional observations. Further- 
more, the improvement in the figure of merit induced by weak 
lensing mass calibration is almost independent of the magni- 
tude of this scatter, and should thus apply to a wide range of 
mass proxies. 
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APPENDIX 
SOURCE CLUSTERING 

In the main section of this paper, we assumed source galaxies are randomly distributed in the sky. We now investigate whether 
source clustering can have an impact on our results. Source clustering implies that equation[lO]must be replaced by 

<^> = ^^4q + ^ (A1) 

where £y is the projected galaxy-galaxy correlation function. Note that in general this term will correlate different radial bins, so 
that if A£ Q and AE^ are the estimators at radial bins a and j3, the new term in the covariance matrix is given by 

C a p = t 2 c ^ 2 nfg agf3 V aP (A2) 

where 

V «P = tV E A! ^ l W ■ (A3) 

' IJ 

and A a and Ap are the area of the annulus a and j3 respectively. Note we have used the thin annulus approximation to set (x and 
g constant within each annulus. Taking now the continuum limit, we arrive at 



y ^ = rV / P(l)\K0Wp(l)\. (A4) 

A a Ap J (Z7r) z 



Before going any further, it is worth taking a second to compare equation |A2 1 with equation [12] Specifically, note that shape 
noise is explicitly dependent on the source density, whereas source clustering is not. This implies that at a sufficiently high source 
density, source clustering must dominate. On the other hand, source clustering scales as g 2 , whereas shape noise scales as a 2 . 
More precisely, source clustering will be relevant if gly aa is comparable to <j 2 ,/2hA. Thus, unless the source density is quite 
large, we expect shape noise to dominate. What follows is a quantitative confirmation of this expectation. 

To do so, we need to begin by estimating V a p, which in turn requires that we compute the projected source galaxy power 
spectrum. Let then n co {x) be the comoving galaxy density field. Assuming flatness, the corresponding projected galaxy density 
field is 

n(ff) = Jdzn coX 2 ^H(z-z L ) (AS) 

where H is a step function that selects only galaxies at redshift larger than the lens redshift of interest zl- The fluctuations in the 
source density field are therefore given by 

5(9) = J d X g(x)6 g (x) (A6) 

where 5 g (x) is the 3D galaxy density contrast, and 

g(x)= rj T , Itt, ;f(d(zr) H(z-zl). (AT) 



Jdzf(z)H(z-z L ) J \dz 

Using Limber's approximation, the angular power spectrum of the source density field is related to the three dimensional galaxy 
power spectrum via 

P(l) = J ' d X ^-P iD (l/x,z). (A8) 

All that remains is to specify the 3D galaxy power spectrum. To do so, we ass ume a constant scale -independent bias b = 1 relative 
to the no n-linear matter po wer spectrum P„, m (k,z), which we compute using S mith et al.1 (12003b as implemented in the CAMB 
package (Lewis et al.ll2000l) . 

Having determined the projected source power spectrum, we need to compute the mean power over the filter functions W a . We 
assume W a is a top hat in radius going from log 10 /? £ [log l0 R c - A,log 10 /? f + A], with A = 0.02 as per our fiducial assumptions. 
The Fourier transform of W a is therefore 

W a (l) = J d 2 e W a (6) exp(;7 -9) = 2ttJ """ dO 9J (W) (A9) 
_2tt 

— ~j2~ \MmaxJ\{Wmax^~ Wmi n J\{u m i n y\ (A.10) 
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FIG. 9. — Ratio of the shape noise error considered in the main body of the manuscript to our estimated source clustering error using our fiducial assumptions 

for the estimator AE. The ratio is shown for halos of mass M = 10 14 Mq (thin lines) and M = 10 15 Mq. Different lines correspond to different redshifts, namely 
Z = 0.1 (solid), z = 0.5 (dashed), and z = 0.9 (dotted). In all cases, we assume a concentration c = 5. We conclude that for the source densities expected for the 
DES source clustering is a subdominant source of noise. 

where u = Rl/d^ and d& is the angular diameter distance to the redshift of the lens. We now have all the ingredients necessary for 
computing the source clustering error. 

Figure|9]compares the ratio between the shape noise term used in the main section of the paper, and the diagonal contribution 
C aa to the covariance matrix due to source clustering. Note we plot the ratio of the errors (i.e. the square root of the variance) 
rather than the ratio of the variance. The ratio of the noise terms is estimated for halos of mass M = 10 14 Mq (thin lines) and 
M = 10 15 M Q (thick lines) for three different redshifts, z = 0.1 (solid), z = 0.5 (dashed), and z = 0.9 (dotted). We can see that 
source clustering is always significantly smaller than the corresponding shape noise terms, so that it can be safely neglected in 
our forecast. It is worth remarking, however, that the ratio of these two noise terms is dependent on the assumed radial bin 
width, with broader bins leading to smaller ratios reflecting the decreased shape noise. Figure [9] assumes our fiducial bin width 
±Alog 10 = 0.02. If we were to use relative broad radial bins of width ±Alog 10 /? = 0.1, the error ratio can drop to a factor of 3 
for the 10 15 Mq halos and down to a factor of ~ 10 for the 10 14 Mq halos. Thus, even when employing broad radial bins, the 
errors are dominated by shape noise rather than by source clustering. 

THE MEAN AND VARIANCE OF As' 

The estimator A S takes the form A E = E e ;c/y where 

x=^AQ(l + 6i)eiWi (Bl) 

i 

y = ^An(l + SiWi. (B2) 

i 

We write x = x+ Ax andy = y+ Ay, and use the binomial expansion to solve for AE assuming Ay/y <C 1 to expand up to second 
order. We find 



AE' = E c - 

y 



Ax Ay AxAy Ay 1 

1 + — f 1 =^- + ^n 

x y xy y 



(B3) 



Upon taking the expectation value, and using the fact that x = Ag and y = A, we arrive at / AE \ = E c g, as per equation 
Squaring equation |B 3 1 we find 

Ax 2 3Ay 2 4 AxAy 



(AE Y = E^ 



1 + ^^ + - 



jr y xy 
To compute the expectation value, we use the fact that x = gA, y = A, and 



(B4) 



W a l* + &)jjk (B5) 
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(AxAy) = 
<A/) = 
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gA 

p q l 1 h 
A 



The expectation value of ( AS ) simplifies to 



(A£) 2 ) = £ 



^2 2 



1 + 



1 < 

[itl 2 nA 2g 2 



(B6) 
(B7) 

(B8) 



which leads directly to equation^ 



A MODEL FOR CLUSTER GALAXIES 



In order to compute this systematic error due to <5r c i we must first estimate the ratio n c \/n. For this, we derive an order-of- 
magnitude estimate as follows: let e be the fraction of galaxies within a narrow redshift slice ±Az about the redshift of the 
cluster, and ±L be the corresponding physical width of the slice. These two quantities are related to the redshift width of the slice 
via 

Idh f(z) 
e= -~r2Az= ^— , 2Az 



ndz f™dz'f(z') 



where f(z) is given by equation[14] and 
Now, the mean three dimensional galaxy density p g within the redshift slice is 



l+z 



2LDl 



Letting S g be the galaxy density contrast field, the projected cluster galaxy density is given by 



dx 



dx 



en f dx 



Assuming galaxies trace mass, we can set 6 m = 6 g , and therefore 

(l + S g ) = (l + 6 m )-- 

Replacing in our expression for n c \, 

en f dx p n 



Pm(z) ' 



2L J \+zpm{zL) 2Lp m (z) 
Inserting our expressions for e and L, and setting p m = il m p c (l +z) 3 we arrive finally at 

«cl _ f(z) E 

n ~ /J 50 dz! f(z') ( 1 + z) 2 n mPc cH-^ (z) ' 



(CI) 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 



(C7) 



